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We consider in detail the nonlinear equations encountered at each time 
step when certain implicit numerical-integration algorithms are used. 
In terms of only the properties oj the Jacobian matrix of the pertinent set of 
differential equations, we present necessary and sufficient conditions for 
the existence and uniqueness of the solution of the nonlinear equations for 
all continuous forcing functions and any given step size. Since engineers 
often think about dynamic nonlinear transistor network problems in terms 
of the eigenvalues of the relevant Jacobian matrix, the results described are 
of immediate conceptual value. In particular, it is possible to carry out 
the algorithms whenever the conditions presented are satisfied. 

Several other types of results are also presented. For example, for a 
special but significant and useful numerical-integration formula, theorems 
are proved concerning properties of the computed sequence such as the 
extent to which the sequence is relatively immune to small local errors in- 
troduced at each step as a result of the fact that it is ordinarily not possible 
to compute the solution of a certain equation exactly. 

All of the results are concerned ivith network models that are often used, 
in computer simulations. In fact, ice heavily exploit some special properties 
possessed by the nonlinear functions associated with such models. 

I. INTRODUCTION 

The set P of all real square matrices each with all principal minors 
nonnegative plays a key role in the study 1-3 of nonlinear equations of 
the form F(x) + Ax = B, and more generally 4 of equations of the form 
CF{x) + Ax = B, in which F{-) is a ''diagonal monotone-nondecreasing 
mapping" of real Euclidean n-space E n into itself, A and C are real 
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n X n matrices and B is an element of E" . Such equations arise in the 
dc analysis of transistor networks, the computation of the transient 
response of transistor networks, and the numerical solution of certain 
nonlinear partial-differential equations. 

In Ref. 3 a nonuniqueness theorem is proved which focuses attention 
on a simple special property of transistor-type nonlinearities. It shows 
that for any transistor-type exponential F(-) the equation F(x) -+- 
Ax = B has at least two solutions x for some B zE 1 whenever A 4 P . 
The theorem shows that some earlier conditions 1,2 for the existence of a 
unique solution cannot be improved by taking into account more in- 
formation concerning the nonlinearities, and therefore makes more 
clear that the set of matrices P plays a basic role in the theory of 
nonlinear transistor networks. Ref. 3 also contains material concerned 
with the convergence of algorithms for computing the solution of 
/P(a>) -{-Ax = B as well as of more general equations, and some related 
problems concerning the numerical integration of the ordinary dif- 
ferential equations which govern the transient response of nonlinear 
transistor networks are considered briefly. 

The primary purpose of this paper is to present the results of a 
continuation of the numerical integration study initiated in Ref. 3. 
Here we further exploit the special property of transistor-type exponen- 
tial nonlinearities used in Ref. 3. 

We consider in detail the nonlinear equations encountered at each 
time step when certain implicit numerical-integration algorithms are 
used, and, in terms of only the properties of the Jacobian matrix of the 
pertinent set of differential equations, we present necessary and suffi- 
cient conditions for the existence and uniqueness of the solution of the 
nonlinear equations for all continuous forcing functions and any given 
step size. Since engineers often think about dynamic nonlinear transistor 
network problems in terms of the location of the eigenvalues of the 
relevant Jacobian matrix, the results described in Section 2.2 are of 
immediate conceptual value. In particular, these results are of a very 
different character than those that appear in the literature, and when- 
ever the conditions presented are satisfied, it is possible to carry out 
the algorithms. Under the assumption that the conditions are satisfied, 
we also show that there are convergent algorithms for solving the non- 
linear equations, and that the Jacobian matrix of the nonlinear equa- 
tions is essentially always at least weakly well-conditioned in a sig- 
nificant sense. 

A part of Section 2.3 reports on a general result concerning conditions 
under which it is possible to invert nonlinear mappings in E" . More 
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explicitly, we show that a proposition proved by G. H. Meyer enables 
us to give a short proof of a new theorem which is a considerably 
stronger result than that described and used in Ref. 11. 

We also present a set of results concerning properties of an important 
class of transistor-diode networks for which certain implicit numerical- 
integration algorithms can be carried out for all values of the step size, 
and, for a special but significant and useful numerical-integration for- 
mula, theorems are proved concerning some properties of the computed 
sequence such as the extent to which the sequence is relatively immune 
to small local errors introduced at each step as a result of the fact that 
it is ordinarily not possible to compute the solution of a certain equation 
exactly. 

Finally, in addition to other results, we present new theorems con- 
cerning the existence of solutions of the nonlinear dc equation under 
very realistic assumptions from the viewpoint of models often used 
in computer simulations. 1 " 

Section II contains a detailed discussion of the results and their 
significance. 

II. TRANSIENT RESPONSE OF TRANSISTOR-DIODE NETWORKS AND IM- 
PLICIT NUMERICAL-INTEGRATION FORMULAS 

2.1 Introduction 

We shall consider explicitly only networks containing transistors, 
diodes, and resistors. However, the material to be presented can be 
extended to take into account other types of elements as well. In addi- 
tion, we shall focus attention on the use of linear multipoint integration 
formulas of closed (i.e., of implicit) type, since such formulas are of 
considerable use in connection with the typically "stiff systems" of 
differential equations encountered. 

A very large class of networks containing resistors, transistors, and 
diodes modeled in a standard manner is governed by the equation 5 "* 

j t + TF[C~\u)) + GC~\u) = B(t), / £ (1) 



t Results concerning the dc equation lire directly relevant to the problem of 
computing the transient response to the extent that in order to numerically integrate 
the differential equations it is ordinarily necessary to first solve a dc problem to 
determine the initial conditions. 

* As a practical matter, the models of transistors and diodes employed here are 
often used in computer simulations. Of course in some cases it is necessary to use 
more complicated models. 
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with G = G(I + R<2)~ 1 and where, assuming that there are q diodes 
and p transistors, 
(i) T = T 1 T 2 • • • T p /„ , the direct sum of the identity 
matrix of order q and p 2 X 2 matrices T k in which 
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with < a { T k) < 1 and < a ( f k) < 1 for k = 1, 2, • • • , p. 
(ii) R = R t # 2 • • • R„ #o , the direct sum of a diagonal 
matrix # = diag (■/■, , r 2 , ■ ■ ■ , r B ) with r k ^ for k = 1, 2, • • • , q 
and p 2 X 2 matrices i?* in which for all A: == 1, 2, • • • , p 



Rl = 



ri k) 



I b 



(I.) 



ri w + r{ w J 

with r< fc) ^ 0, r[ h) ^ 0, and rj ft) ^ 0. (The matrix 72 takes into 
account the presence of bulk resistance in series with the diodes 
and the emitter, base, and collector leads of the transistors.) 
(Hi) (t is the short-circuit conductance matrix associated with the re- 
sistors of the network. (It does not take into account the bulk 
resistances of the semiconductor devices.) 
(iv) F(-) is a mapping of E lZv * Q) into # (2p+0) defined by the condition 
that 

F(x) = [fiM, U(x 2 ), ■■■ , / 2p+(J (z2 P+0 )] tr 

for all x e E (2v+Q) with each /,(•) a continuously-differentiable 
mapping of E l into E 1 such that /;(«) > for all a t E l . 
(«) C'\-) is the inverse of the mapping C(-), of # (2p+0) into itself, 
defined by 

C(x) = ex + tF(x) 

for all x t E i2p+a) with c = diag (c, , c 2 , • • • , c (2p+ „), t = diag (tj , 

T2 , • • • , 7" ( 2„ + <,)), and with each t, and each c, a positive constant. 

(yi) B(t) is a (2p + g)-vector which takes into account the voltage and 

current generators present in the network, and 
(dm) u is related to v the vector of ideal-junction voltages of the semi- 
conductor devices (v does not take in account the voltage drops 
across the bulk resistors) through C(v) = w for all v t E {2j " . 

Equation (1) is equivalent to + 

t In Ref. 5 it is shown if B(- ) is a continuous mapping of [0, °0 into 2j?<»p + <»>, then 
for any initial condition u (0 > z E^^ there exists a unique continuous (2p + q)- 
vector-valued function u(-) such that u(0) = u<°» and (1) is satisfied for all i > 0. 
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u + f(u, t) = B i2p+Q) t ^ (2) 

in which of course 

f(u, t) = TF[C-\u)\ + GC-\u) - Bit) (3) 

and 0<2,> k,) is the zero vector of order (2p + q). 

It is well known that certain specializations of the general multipoint 
formula 8,7 

r r 

.'/„.! = £ O-klJn-k + h J2 b k yn-t ( 4 ) 

fc-0 A=-l 

in which 

//„-, = -/(?/„-* , (n - k)h) (5) 

can be used as a basis for computing the solution of equation (2). Here 
h, a positive number, is the step size, the a k and the b k are real numbers, 
and of course y n is the approximation to u(nh) for n ^ 1. 

In the literature dealing with formulas of the t3'pe (4) in connection 
with systems of equations of the type (2), information concerning the 
location of the eigenvalues of the Jacobian matrix J u of j{u, t) with re- 
spect to u plays an important role in determining whether or not a given 
formula will be (in some suitable sense) stable. In particular, an as- 
sumption often made is that all of the eigenvalues of J u he in the strict 
right-half plane for all t ^ and all u. For f(u, t) given by equation (3), 
we have 

./„ = T diag ( tffrffi f)1 , J + d diag { , ),. , r~\ (6) 

in which for j = 1,2, • • • , (2p + g) (7,(m,) is the jth component of C~ l (u). 
Thus here J u is a matrix of the form 

TD t + GD 2 (7) 

where Z>, and D 2 are diagonal matrices with positive diagonal elements. 
A simple result concerning (7), Theorem 4 of Ref. 3, asserts that if there 
exists a diagonal matrix D with positive diagonal elements such that 1 

(»") DT is strongly column-sum dominant, and 
(ii) DO is weakly column-sum dominant, 
then for all diagonal matrices D, and D 2 with positive diagonal elements, 

f The terms "strongly-eolumii sum dominant" and "weakly column-sum domi- 
nant" are reasonably standard. However, I hey are defined in Section III. 
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all eigenvalues of (7) lie in the strict right-half plane. This condition 
on T and G is often satisfied.* 

The subclass of numerical integration formulas (4) defined by the 
condition that b_j > are of considerable use 8 " 10 in applications involv- 
ing the typically "stiff systems" of differential equations encountered 
in the analysis of nonlinear transistor networks. With 6_ t > 0, y n +j is 
defined implicitly through 

r r 

y n+1 + hb^i(y n+1 , (n + l)h) = £ a k y n - k +ftl b k y n . k 

*=o *- o 

in which the right side depends on y n - k only for k e {0, 1, 2, • • • , r}, and 
for f(u, I) given by equation (3), we have 

y n+l + Kb^lTFiC- 1 ^)] + OC' l (y n+l )\ = q„ (8) 

in which 

k=0 k=0 

Obviously, the numerical integration formula (8) makes sense only if 
there exists for each n a y n+1 t E iZp+9) such that (8) is satisfied. 

2.2 The Jacobian Matrix J u and Necessary and Sufficient Conditions for 
the Existence of a Unique Solution y n+1 of (8) for All q n tE 2l 
Here we shall make the additional assumption that the functions 
/,(•) are such that the mapping F(-) belongs to the set JFj ap+a) defined 
in Section 3.1. This assumption is satisfied whenever the /,(•) are the 
usual Ebers-Moll exponential-type nonlinearities. That is, $ l 2p+Q) con- 
tains all of the mappings F( • ) such that for each •/' 

f s {Xi) = a,[exp (b/i,-) - 1] or f t (x t ) = a,[l - exp (-b t x t )] 

for all Xj t E 1 with a,- and 6, positive constants. 

Our first result, Theorem 1 of Section III, is a rather strong result 
concerning the relation between properties of the Jacobian matrix J u 
and properties of equation (8). Let 3 denote the set of all real numbers 
a such that det (<xl + J u ) = for some u t E <2p+Q) . In other words, let 
S denote the set of all real numbers a such that —a is an eigenvalue of 
J u at some point u. According to Theorem 1, equation (8) possesses a 
unique solution y a+1 for each q n t E l2p+Q) (and hence eachJ5[(n + l)/t] e 
E i2v+Q) ) if and only if (Aft-i)"" 1 t S, and also if (M-i)" 1 1 S then equation 
(8) possesses at least two solutions for some q n t E vlp+Q) (and hence for 

f See Ref. 5 for examples. 
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some B[(n + l)/t] t E lap * 9) ). Therefore, in particular, equation (8) pos- 
sesses a unique solution for all q n t E l2p+9) and all h e (0, h], in which h 
is an arbitrary positive constant, if and only if the intersection of the 
interval [(ftft-i)" 1 , <*>) and S is the null set, and equation (8) possesses 
a unique solution for all q n t E {2p *" ) and all h > if and only if S contains 
no points of the interval (0, «). Finally, as a somewhat peripheral 
matter, according to Theorem 1, the dc equation TF{v) -\- Gv = B has 
at most one solution v for each B i E (2p+a) if and only if i %. 

The statements made in the preceding paragraph are surprising to 
the extent that on the one hand they are rather definitive and on the 
other hand they involve only the location of the real eigenvalues of 
«/"„ .* Since engineers often find it helpful to think about nonlinear 
systems in terms of the location of the eigenvalues of a pertinent 
Jacobian matrix, it is also of interest to note here that equation (8) can 
possess more than one solution y n+l for some q a and some h > only if 
the transistor-diode network is locally exponentially unstable at some 
operating point, that is, only if at some operating point u, —J u has a 
real positive eigenvalue. 

2.3 Existence of Convergent A Igorithms for Computing the Solution of (8) 

Throughout this section we assume that the /,(•) are such that the 
additional condition that F(-) t 5 ( 2p+Q) is satisfied. 

Whenever (7i6_,) -1 is not contained in the set S of Section 2.2, equation 
(8), which we shall write as Q(y n +i) = q n , possesses a unique solution 
y n+i for any q n t E i2v+Q) . We show here that when (/i&_,) _1 ^ S and each 
/,-(•) is twice continuously differentiable on E 1 ,* there exist steepest 
descent as well as Newton-type algorithms each of which generates a 
sequence in E {2l,+Q) which converges to y n - x . 

Assume that (hb-x) i 2. The Jacobian matrix (/ + hb-xJ^^) of Q(-) 
satisfies 

det (/ + /ife-,J r/n+1 ) ^ for all y n+l e E l2p+Q) . (9) 

Hence Q(-) is a local homeomorphism on E V2p * q) and since there exists 
a unique y n+1 £ E (2p + Q) such that Q(y n+l ) = q n for each q n z E Uv+,,) , Q(-) 



t Indeed, while we can write (8) as Q(y n+ i) = q n with Q(-) a conlhiuously-dif- 
ferentiable mapping of Zs'< 2 ' ,+ fl> into itself with Jacobian matrix (7 + /ife_,./,, n+i ) recall 
that for R(-) a general continiiously-differentiable mapping of E n into itself with 
Jacobian matrix J, det .7^0 throughout E n does not imply that (and is not implied 
by the statement that) for each x e E" there exists a unique y e E n such that 
R(y) = x, even for n = 1. 

J This differentiability condition is obviously satisfied if the /,(•) are the usual 
exponential functions. 
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is a homeomorphism of E (2v+a) onto itself. Thus, with ||-|| any norm 

«~ ip(2p + a) 

on cj , 

II Q(y) II — > °° as H y II -» w. f 

Let #(•) be defined by the condition that R(y) = Q(y) — q n for all 
y c E i2v * q) . Then fl(-) satisfies || R(y) || -> «> as || ?/ || -» « and the 
determinant of the Jacobian matrix of R(-) does not vanish throughout 
2£<2p+a> ^ Therefore, assuming that R(-) is twice continuously differenti- 
able on E (2p+Q) , it follows (see the Appendix) that the solution y n+1 of 
R(y n+1 ) = (2p+a ) can be computed by using certain steepest descent or 
Newton-type algorithms. 

2.4 The Jacobian Matrix (I + /i6_iJ„ n+1 ), and Inversion oj Nonlinear 
Operators on W and Jacobian Matrices 

As in Section 2.3, let the additional condition that F(-) c ff£ 2p+ * ) be 

satisfied and let Q(- ) be the mapping of # (2p+0) into itself with the prop- 
erty that equation (8) can be written as Q(y n+1 ) = q„ . According to 
Theorem 2 of Section III the Jacobian matrix (/ + hb-iJ Un+i ) possesses 
the property that there exists a constant e > such that 

det (1 + to-iJ^) £ e for all y n+1 1 E tap+9) (10) 

if and only if the matrix 

[(M-iTV + rrw-o-'o + g]i 

which we shall call *S, belongs to the set P of all real square matrices 
each with all principal minors positive. Thus when S e P the matrix 
(Z + hb-iJ Vn+1 ) is well conditioned in at least the weak sense of (10). 
This fact is of some interest for two reasons. First, certain standard 
algorithms require that the matrix (I + ft&-i«7» ll+1 ) be inverted along a 
sequence of points \yl+\ } in order to compute the solution y n+i of equa- 
tion (8), and, secondly, Theorem 3 of Section III shows that if 
det [(/i6_,)"7 + Jul * for aU u c E V2p+tl) and all (M-,)" 1 t d' in which 
d' denotes either (0, «>) or any interval contained in (0, «>), then S c P 
for all but at most a finite number of points (hb.i)~ l contained in S'. 
Therefore, referring to the material of Section 2.2, if Q(y n +i) = q n pos- 
sesses a unique solution y n+1 for all q n t E <2j,+Q) and all (/ib_i) _1 e d', then 
(I + hb-iJ yu+l ) is at least weakly well conditioned at all but at most a 
finite number of points contained in 6' . 



f Since Q( • ) is a homeomorphism of j£< 2 » + «> onto itself, Q{ ■ ) -1 exists and is contin- 
uous. Therefore, the image of any closed ball in Jjto+e) under Q(- ) _1 is contained in 
some closed ball in E (2 r +9) , and hence || Q(y) || -> » as \\y || -» =>. 
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Since the elements of (/ -f- hb^iJ Vn+x ) are bounded on y n+1 e 2£ (ap+a) j ft 
follows from a theorem described by M. Vehovec 11 that for each q n t 
E l2p+a) there exists a unique y n+1 e E l2p+a) such that Q(y n+1 ) = g n if S e P. 
More explicitly, the theorem described 1 by Vohovec asserts that if R(-) is 
a continuously-differentiable mapping of E 1 into E" with J(R) a the Jaco- 
bian matrix of R(-) at an arbitrary point q e E 1 , if the elements of 
J(R)„ are bounded on E n , and if there exists a positive constant e such 
that det J{R) q ^ e for all q t E", then R(-) is a homeomorphism. Thus, 
using the theorem of Ref. 11 and Theorems 2 and 3 of Section III, 
we are able to show that if det [(A6-0" 1 * + /„] ^ for all u z E i2v+q) 
and all (hb-,)~ l t 6' , then for all but at most a finite number of points 
(A6_i) -1 t H', (8) possesses a unique solution y n+1 for each q n e £ ,(2p+<,) . 
Although this result is obviously much weaker than the existence 
proposition presented in Section 2.2, it shows that the theorem of 
Ref. 11 can be exploited to provide some insight in connection with 
the specific problem considered here. 

The theorem of Ref. 11 is of interest primarily because the key hy- 
pothesis concerns only the determinant det J(R) a (as opposed to the 
condition of Palais 1 that || R(q) || — > oo as \\ q\\ — * °°). Theorem 4 of 
Section III is a general result which is considerably stronger than the 
theorem of Ref. 11. It shows that the condition of the theorem of Ref. 
11 that there exist a positive constant e such that det J(R) a ^ e for all 
q can be replaced with the condition that there exist real constants a > 
and 6^0 such that 

det J{R) a > — —TTx — IT f o r all q t E". 
a + o I) q (I 

2.5 A Class of Netiuorks for Which (8) Possesses a Unique Solution for 
All Values of the Step Size 

There is an interesting class of transistor-diode-resistor networks 
with the property that for each network in the class, equation (8) pos- 
sesses a unique solution for all h > (i.e., for all h > 0, all q n t E i2v+Q) , 
and all diagonal matrices c and t with positive diagonal elements). In 
order to define and discuss that class, consider the dc equation TF(v) -f- 
Gv = B in which v is the (2p + q) -vector of semiconductor ideal-junction 
voltages and B t E l2p+a) . If p > and the matrix R of Section 2.1 is 
the zero matrix, v, is the emitter-to-base voltage of transistor one, v 2 
is the collector-to-base voltage of transistor one, and so forth. By port 

1 According to Vehovec, the theorem was recently proved by I. Vidar, and the 
proof is expected to appear in the journal Glasnik Matematicki. 

* See Ref. 12 and the appendix of Ref. 13. Here ||-|| denotes any norm on E". 
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j of the transistor-diode-resistor network we mean the terminal pair 
between which the voltage v, appears. Again we shall make the assump- 
tion that F(-) t ^ 2p+a) . 

In Ref. 3 it is proved that TF(v) + Gv = B possesses at most one 
solution v for each B t E l2p+Q) if and only if T~*G c P„ . It is also proved 
in Ref. 3 that equation (8) possesses a unique solution y n+1 for each q n t 
E (2v+a) and each h > if M~ l G z P for all M z 3(T) in which here 3(T) 
denotes the set of all real matrices having the same form as T and with 
the "a's" of M not larger than those of T. f In other words, it was also 
proved in Ref. 3 that equation (8) possesses a unique solution ?/ n+1 for 
each q n z E i2p+Q) and each ft > if the dc equation possesses at most 
one solution for each B z 2? (2p+< ' ) for "the original set of a's as well as 
for an arbitrary set of not-larger a's." Before proceeding, and for the 
sake of completeness, we mention here that the same result can be ob- 
tained by way of the approach of Section 2.2; a direct corollary of The- 
orem 5 of Section III, Corollary 1, shows that if M" l G z P for all M z 
3(9"), then det (trl + J u ) ^ for all real a ^ and all u z E l2p+Q) . 

Theorem 5 of Section III provides considerable information concerning 
the nature of the class of networks for which M~*G z P for all Me3(T). 
In particular, the theorem shows that M~ X G z P for all M z 3(T) 
if and only if M~ l G e P for all M e 3„(T) in which 3 (T) is the set of 
all 2 2p real square matrices M having the same form as T and with 
each "a" of M either zero or the corresponding "a" of T. The the- 
orem also shows that "M^G t P Q for all M t 3(T)" is equivalent to 
each of six other statements involving T and G. For example, according 
to Theorem 5, we have M~*G e P for all M e 3(T) if and only if either 
r -1 (G -f- D) t P Q for all diagonal matrices D with positive diagonal 
elements, which has an obvious network interpretation in terms of the 
addition of resistors to the network characterized by G, or T~ l G e P Q 
and (TJ^Gy, e P for all pairs of matrices T w and G w obtained from T 
and G, respectively, by deleting an arbitrary set w of rows, and the same 
set of columns, of both T and G. 

When the matrix R of Section 2.1 is the zero matrix, the last condition 
on T and G of the preceding paragraph also has a simple network 
interpretation: Given T and G, we have T~ l G z P , and any network 
obtained from the network characterized by T and G by short-circuiting 
an arbitrary set w of at most all but one of the (2p + q) semiconductor 
junctions possesses the following property. With respect to the voltage 
vector v w associated with the junctions not short-circuited, and with 



♦ See Definition 4 of Section III for a precise definition of §(T). 
1 See Definition 5 of Section III for a precise definition of So{T). 
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the components of v w taken in the same order as those of v, the "new 
T and G" matrices 1 T w and G a satisfy (T^G,, t P . As reasonable 
as this condition or any of the other seven equivalent conditions of 
Theorem 5 might seem, and even though, as Theorem 6 of Section III 
shows, T" l G e P implies that (T„)~ l G w e P whenever w has the property 
that if the port number associated with one junction of a given transistor 
is contained in iv, then the port number associated with the other 
junction of that transistor is also contained in w, it is the case that 
there are transistor-diode-resistor networks for which T~*G e P and 
M~ l G i P for some M t 3(T). In fact, Ref. 14 presents an example in 
which p = 3, q = 0, T~ X G t P , and T~\G + D)tP for some diagonal 
matrix D with positive diagonal elements. However, the class of networks 
for which T~*G t P implies that M~ X G t P for all M t 3(T) is clearly 
quite large; it obviously includes all networks in which p = 0, it includes 
all networks in which the base terminals of all transistors are connected 
to a common point, and as Theorem 7 of Section III shows, the class 
includes all networks in which T~ l G t P and p = 1 or p = 2. n 

2.G Results Concerning the Numerical-Integration Formula y n+i = y n -f- 
hy n -,i 

The general multipoint formula (4) reduces to the well-known 
implicit numerical-integration formula y n+ i = y n + hy n +i when a = 
b-i = 1, b = 0, and a k = b k = for k = 1,2, • • • ,r. For that important 
special case, and with y n+l given by equations (3) and (5), \y n +i] is 
denned implicitly through 

y n+1 + h\TF[C- l {y n ^)] + GC~ l {y n+1 )\ = y n + hB n (11) 

for all n ^ 0, in which B n = B[(n -f- l)h]. Here we describe some detailed 
results concerning the relation between the sequences \y n +i\ and \B n \. 
We assume throughout this section that G is such that there exists a 
diagonal matrix D with positive diagonal elements with the property 
that both DT and DG are strongly column-sum dominant. This condi- 
tion, which is often satisfied, guarantees that there exists a unique 
solution* ?/ n+1 of equation (11) for each (y n + hB n ) t £ ,(2p+a) . 

+ It. is a simple matter to show that the "new T and G" matrices are T w and (?„. 
ft It is proved in Ref. 14 that, if q = and if p = 1 or p = 2, then T~ l G e F 
implies that T~ 1 (G + D) c P for all diagonal matrices with positive diagonal ele- 
ments. Thus, by the equivalence of statements (t) and (v) of Theorem 5 of Section 
III, it follows at, once that if T~ l G c P then M~ l G c P for all M z Q(T) if q = and 
p = 1 or p = 2. The proof of essentially the same end result given here is of a very 
different nature and is quite short. 

5 See Ref. 5 for examples. 

* A result mentioned in Section 2.1 implies that if DT and DG are both strongly 
column-sum dominant, then det {(h)' 1 ! + J u ] ?* for all it c _E< 2 p + 9> and all h > 0. 
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Let || ■ ||, be defined by the condition that || v ||, = £i- 1 + " ) | v, | for all 
v e E i2v * Q) . According to Theorem 8 of Section III, there exists a positive 
constant 8 depending only on the c, , the r,- , T, G, and D such that 

|| Dy n ||, £ (1 + 8hy || Dy [|, + h £ (1 + 5ft)"* || D5 (B _ t , ||, 

for all n ^ 1. Therefore, it follows that for all h > 0, the sequence 
Vi > Ife j • • • is bounded whenever the sequence 5, , J5 2 , • • • is bounded, 
and 7/1,2/2, • • • approaches ( 2 P+9 ) the zero vector of E (2 ' " whenever 
B x , B 2 , • ■ • approaches (2 p+ a ). 

Typically at each step an iterative algorithm is employed to compute 
the solution ?/„+, of equation (11). Since it is ordinarily not possible to 
compute y n +i with infinite precision, it is important to consider the 
effects of the errors which are introduced. While, ideally, we would like 
to determine the sequence \y n +i} defined by equation (11) and some 
initial-condition vector y , suppose that we determine instead a sequence 
{$ n +i } such that, with e an arbitrary positive constant, \\ D(y„ — y*) ||, 
^ e for all n ^ 1 and 

y* +1 + h{TF[C-\y* +l )] + CC"'^* ,)! = y n + hB n (12) 

for all n ^ 0. That is, suppose that at each step the local error || D(y n — 
y*) ||, in solving for "y„ +l " is at most c. Then, according to Theorem 8, 
and with 8 the positive constant referred to above, 

|| D(y n - y n ) ||, ^ (1 + 8hY n \\ D(y - y n ) ||, 

+ « Z (1 + 8hy k for all n ^ 1 

fc-0 

in which y is the approximation to y . Therefore, given an arbitrarily 
small positive constant p, for any h > it is possible to choose y and 
e > such that the accumulated-error vector (y„ — ?/„) satisfies 1 1 y n — 
y n Hi ^ p for all n ^ 1. 

Finally, Theorem 9 of Section III provides us with a conceptually 
interesting uniform bound on the norm of the difference between cor- 
responding elements of the sequences \y n \ and jw„} in which u n = u(nh) 
for all n ^ and m(-) satisfies the differential equation (1). According 
to Theorem 9, there exist positive constants 5 and p, both independent 
of h, such that 

II D(u n - y n ) ||, ^ (1 + 8hY n || D(u Q - 7/0) Hi + ph 

for all n ^ 1, assuming that the elements of B(-) and (d/dt)B{-) are 
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bounded and continuous on [0, »). In particular, if y = u we see that 
there exists a positive constant p, independent of h, such that || u n — 
y n \\\ ^ p'h for all n ^ 1, provided only that the assumptions of this 
section are satisfied and that B{-) and (d/dt)B(-) are bounded and 
continuous on [0, °°). 

2.7 Conditions Which Imply That T~ l G(I + RG)~ l z P 

In this section and in Section 2.8 we present some results concerning 
properties of the dc equation TF(v) + Gv = B. These results are directly 
relevant to the problem of computing the transient response of tran- 
sistor-diode networks to the extent that in order to numerically integrate 
the differential equation (1) it is ordinarily necessary to first solve a 
dc problem to determine the initial conditions. 

As indicated in Section 2.1, G = G(I + RG)' 1 in which R takes into 
account the bulk resistances associated with the semiconductor devices. 
Here we present some material concerning conditions which imply 
that T~ l G{I + RGy 1 belongs to P . 

Let p > 0. Theorem 10 of Section III asserts that T~ X G(I + ^) _1 8? 
whenever T~ l G z P n and R satisfies 

al k) (l - a l r k YY t k) = r ( b k) 

a\ k) (i - a j k) yv e k) = ri k) 

for k = 1, 2, • • • , p. This rather special result shows that if F(-) satisfies 
the additional condition that F(-) belongs to the set 9 r o 2p+( ' ) defined in 
Section 3.1, and if the network associated with T and G possesses the 
property that there is at most one solution v of the dc equation TF(v) + 
Gv = B for each B e J5 (2p+a) ) then it is always possible to add certain 
resistors of positive value in series with each transistor lead such that 
the dc equation of the resulting network possesses at most one solution. 

Theorem 11 of Section III directs attention to the fact that there is 
a nontrivial class of transistor networks for which T~ l G(I + RG)' 1 e P 
for all R. According to Theorem 11, if p > and G is such that T^G z P 
for all "ffl's" (i.e., for all a l r k) and af belonging to (0, 1)), then for any 
particular set of "a's" T~ X 6(J + RG)~ l z P for all 72.* 

Given T, an interesting characterization of the class of short-circuit- 
conductance matrices G such that M~ l G z P for all M z 3(T) is provided 
by Theorem 12 of Section III. 1 According to Theorem 12, M~ l G c P 
for all M z 3(T) if and only if T~ X G{I + RG) " l z P for all R satisfying 
certain inequality-type conditions. In particular, if the base-lead 

f A similar result is proved in Ref. 2 under the assumption that G is not singular, 
* The set $(T) is described in Section 2. . 
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resistance of each transistor is taken to be zero, then M~ l G t P for all 
M e 3(T) implies that T _1 (?(J + RG)~ l t P for all nonnegative values 
of each emitter-lead resistor and each collector-lead resistor. 

2.8 Ebers-Moll Models and the Existence oj a Solution o] TF(v) + Gv = B 

In Section III, a set ff 3 of mappings F(-) is defined such that each 
element of 5 S possesses certain important properties possessed by an 
arbitrary F( ■ ) of the type that arises when an Ebers-Moll exponential- 
nonlinear-function model is used for each transistor and diode. In 
contrast with the set of all F(-) such that each /,(•) is a strictly-mono- 
tone-increasing mapping of E 1 onto E 1 , an arbitrary element F(-) of 
SF 3 possesses the properties that for each ;', /,(•) is bounded on either 
[0, oo ) or (— oo , 0], and the two nonlinear functions associated with the 
same transistor are both bounded on either [0, <») or (— «>, 0]. The set 
5 3 is contained in 5^ 2p+0> and contains every Ebers-Moll exponential- 
nonlinear-function-type F(-). 

The first part of Theorem 13 of Section III asserts that the equation 
TF(v) + Gv = B possesses a unique solution v for each F(-) t5 3 and 
each B t £ (2p+a) if and only if T~'G t P and det G ^ 0. It is the "only if" 
part of this proposition which is the new result presented here. The 
proof exploits some special properties of transformerless resistor net- 
works; it shows that if T~ l G t P but det (7 = 0, then there are functions 
t(-) and d(-), both functions talcing on only the values 1 or —1, such 
that there is no solution v of TF (v) + Gv = B for some B t E {2p+a) for any 
set of Ebers-Moll-modeled transistors and diodes with the property 
that for all k transistor k is a pnp device (as opposed to a npn device) 
if and only if t(k) = 1, and for all j diode j is a p-n junction if and only 
iid(j) = l. + 

The discussion of the preceding paragraph concerning the proof of 
Theorem 13 shows that it is not possible to make stronger assertions 
concerning the existence of a unique solution of TF(v) + Gv = B for 
all B t E l2p+q) for Ebers-Moll-modeled transistors and diodes unless 
we take into account more information about the nature of the semi- 
conductor junctions. A good deal of progress in this direction has 
recently been made, and we state here without proof the following 
complete result dealing with diode-resistor networks. 

Theorem 1J+: X Let p = and q > 0. Let F(-) t S 3 (see Definition 12 o/ 



t In contrast, the proof of the "only if" part of Theorem 3 of Ref. 1 shows that if 
A t P then there is a mapping F( ■ ) with each/y( • ) a linear function such that F(x) + 
Ax = B does not possess a unique solution for all B c E n . 

* The proof of Theorem 14 will be presented in a subsequent paper. 
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Section 3.31), and for j = 1, 2, • ■ ■ , q let s ; equal either 1 or — 1 depending 
on whether f ,(• ) is bounded on [0, ») or (— » , 0], respectively. Then, with A 
any real symmetric nonnegative-definite matrix of order q, there exists a 
unique solution v of F(v) + Av = B for all B e E" if and only if there is no 
real q-vector t) such that jj y^ q , Ai\ = 8 Q , and i\ e S, in which 

S = \y: y £ E" and y iSj ^ for j = 1, 2, • • • , q]* 

III. THEOREMS AND PROOFS 

3.1 Notation and Definitions 

Throughout Section III, 

(i) unless stated otherwise, p and q denote nonnegative integers 
such that (p + q) > 0, and n denotes an arbitrary positive 
integer; 
(ii) the set of all real n-vectors is denoted by E", 6 is the zero element 
of E n , and if v t E" and j is an integer such that 1 ^ j ^ n, then 
Vj denotes the jth component of v; 
{Hi) || v || = (]C"=i v ?) 1/2 ail( i II v Hi = 2f-i I v i I for all we 2?" ; for any 
real n X n matrix M, \\ M \\ denotes sup {m: \\ Mx \\ ^ m \\ x \\, 
xtE"}; 
(iv) the transpose of an arbitrary (not necessarily square) matrix 

M is denoted by M tr ; 
(v) I n denotes the identity matrix of order n, and / denotes the 
identity matrix of order determined by the context in which 
the symbol is used; if Qi , Q 2 , • • ■ , Q n are square matrices, then 
Qi Q2 © • • • © Qn denotes the direct sum of Q, , Q 2 , • ■ ■ ,Q n , 
in the order indicated; 
(vi) if D is a real diagonal matrix, then D > 0(D ^ 0) means that 
the diagonal elements of D are positive (nonnegative); and 
(vii) we say that a real n X n matrix M is strongly (weakly) column- 
sum dominant if and only if for j = 1, 2, ■ ■ ■ , n 

m,i > (^) 2 I nti, |. 

Definition 1: The set of all real square matrices il/ such that every princi- 
pal minor of M is nonnegative (positive) is denoted by P (P)- 

Definition 2: Let 3 : o 2p4 ^ ) denote that collection of mappings of E Uv+Q) 
into itself defined by: F(-) t < 5 i 2p+,1) if and only if there exist for ; = 



J In the network case, A = G, and it is often possible to determine by inspection 
whether or not there exists an rj 9^ Q such that Grj = 0, and 77 e S. 
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1, 2, ••• , (2p + q) continuous functions /,(•) mapping E 1 into E l 
such that for each x z # (2p+c) , F(x) = [/, fa), Ufa), • • • , / (2j>+ „ (2p+a ,)] tr , 
and 

(i) inf [/,(« + /3) - /,(«)] = 0, 

ot(-oo.oo) 

(it) sup [/,-(« + /3) -/,-(«)] = +» 

ac(-W.oo) 

for all > and all j = 1, 2, ■ • ■ , (2p + ?)• 

Definition 3: Let 3 denote the set of all real matrices M such that M = 
M, © M a © • • • © M p © /„ with 



M, = 



1 -aj 

( 

L— «/ 



(A')" 



.(*) 1 



^ ai" < 1, and ^ a}" < 1 for all k = 1, 2, • • • , p. As suggested, 

if g = 0, then M = M x © M t © • • • © M p , while if p = 0, then M = 

/.. 

Assumption 1: Throughout Section III, (7 denotes a real nonnegative- 

definite matrix of order (2p + q) . 

A tool that we shall use often is: 

Lemma 1: A real square matrix M is an element of P if and only if det 
(D + M ) 5^ for all real diagonal matrices D > 0. 

Lemma 1 is proved in Ref. 2. 

3.2 Theorem 1: Let F(-) t 5 l 2p * a) with each /,(•) continuously differentiable 
on (- oo, oo ) and #(«) > for all a t (-», »). Lei r e 3, fet C(-) 
[tfiaZ is, c + t^(-)L G, and J u be as defined in Section 2.1, and let a be a 
real nonnegative constant. Then 

ay + TF[C-\y)] + GC'^y) = r (13) 

possesses at most one solution y for each r t E (2p+Q) if and only if 

det (al 4- J u ) ^ for all u t E V2p+9) , (14) 

and if a > and condition (14) is satisfied then for each r t # (2p+<,) there 
exists a solution y of (13). 

3.3 Proof of Theorem I 
We have 
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det (al + ./„) 

= det (*I + TF'[g(u)]{c + rF'[g(u)]}- 1 + G{c + r^'^)]}- 1 ) 

= det \c + rF'[g(u)]}- 1 -det [ac + crF'[g(u)) + TF'[g(u)) + (?} , 

in which g(-) is the mapping of E (2p+Q) onto itself defined by g(u) = 
C-\u) for all u z # (2p+0 \ and F'[g(u)] = diag [/,%,(«,)]}. Since det 
|c + tF'[o(w)]) > for all u, det (<r7 + /„) ^ for all w if and only if 

det{((TT + T)F'[g(u)} + (<rc + G)\ ^ for all t*. 

For each ;' g t (-) maps E 1 onto i? 1 , and since F{-) z 5 ( 2p+a) with each 
/,(•) continuously differentiable on (— co, co) and /,-(a) > for all 
a z (— oo, co), the image of E l under the mapping /,%,(•)] is (0, °°) + 
for all j. Thus, by Lemma 1 (since det(ar + T) 9* 0)(<tt + T)~ l {ac + 
(?) t P if and only if 

det(a7 + ./„) ^ for all u. (15) 

The equation 

ay + TF[C- X {y)] + GC~ l (y) = r 

possesses a solution y if and only if x = C~ l (y) satisfies 

<tC{x) + TF(x) + Gx = r, 

that is, if and only if 

(ar + T)F(x) + (ac + G)x = r. (16) 

But equation (16) possesses at most one solution for each r z E i2v+a) 
if and only if (ar + T)~ l (pc + 6*) c P (see pp. 10.5-107 of Ref. 3) and 
hence if and only if condition (15) is met. 

Suppose now that a > 0. Since G is nonnegative definite, det(<7C + 
(?) * 0. If condition (15) is satisfied then (ar + T)~\ac + (?) z P 
and hence for each r z E l2p+Q) , equation (16) possesses a solution x 
(seep. 99 of Ref. 3).D 

3.4 Theorem 2: Let T z 3, and Ze< F(-) z 5 ( 2p+q) with each /,(•) continuously 
differentiable on (— oo, co) a/»d /J (a) > jor all a z (— co, co). 77ie?i /or 
each a ^ Mere eziste a positive constant e s?*c/j Ma/ det(<rJ + •/.,) ^ e 
/or aZZ w z E i2 " +Q) if and only if (ar + T)~\ac + (?) c P. 



t For any > and any a c (— », » ), //(a + /3) — /y(o) = /3//(5) for some 
5 t [a, a + 0]. 



1756 THE BELL SYSTEM TECHNICAL JOURNAL, OCTOBER 1070 

3.5 Proof of Theorem 2 
We have 
det (<rl + J„) 

= det {d + TF'[g(u)]{c + rF'Mu)}}-' + G{c + tF'^u)}]- 1 ) 
= det \c + rF'[ff(w)]} _1 'det {{or + !T)*"[ff(u)] + (^ + G)\ 

= detiaT + T) j^mi±AL. (17) 

II fc, + r,«toi(«i)D 

in which A = (<rr + r) -1 (i« + G). 

For each sequence e l , e 2 , ■ •■ , e (2p+e ) with each e, either zero or unity 
and e x , e 2 , ■ ■ ■ , e (2 p +a ) not the sequence 1, 1, • • • , 1: let m,,,,,, .... 

«(.p + «) denote the determinant obtained from A by deleting rows 
Pi } f>2 , • • • t Pi an d columns pi , p 2 , • • • , P/ in which ! pi , p 2 , • • ■ > Pj I = 
{;: e,- = l\. Thus for each sequence e x , e 2 , • • • , e (2p+a) other than the 
sequence 1,1, • ■ • ,1 n*, t . «......• <„+„ is a principal minor of A. Let 

TOi.i....,, = 1, and let d t = /,%,(«,)] for all j. Then by a standard 

expression 15 for the determinant of the sum of two matrices 

det (F'[g(u)] + A) = £' d'S<%' ••• cT^Jiii.,.......,.,.^, 

in which ^' denotes a summation over all 2 ( " p+< ' ) sequences e Y , e 2 , • • • , 
e (2p+1J ) and d° = 1 for all j. It is clear that 

j-i 
in which each c,,,,,,...,, ,,,+„ is a positive constant. Thus with 77 = 
det(ar + T), 

Suppose that all principal minors of A are positive. Then there is a 
positive constant 5 such that 

m.,..,.-,. ( „ t ,i = &:».>•. — .M»+f) 
for all e, , e 2 , ■ • • , e (2 p+<,) and hence (since rf, > for all j) det(o7 + 
/«) ^ 7 ? 5forallwe^ (2p+a) . 

As in the proof of Theorem 1, the range of each d t = /-[^(m,)] is (0, »), 
and for any positive constants pi , p 2 , ■ • ■ , p<*»+«) there exists a u t 
E l2p+g) such that d, = p, for all j. If A 4 P then at least one principal 



NONLINEAR NETWORKS 1757 

minor of A is not positive. If A 4 P , then det(F'[g(u) + A}) = 
for some u. Therefore to complete the proof it is sufficient to show that 
if A t P a but A i P then there is no constant e > such that det (a I + 
Ju) > « for all u. 

With A t P and A i P, for at least one sequence e\ , e' 2 , ■ ■ ■ , e' (2p+Q) 

nv,.«',... -..'(.,+.) = 0. 
If det /l = wi .0.....0 = we have 

inf det (a I + J,,) = 

since det(o7 + J u ) — » as (/, — > for all j. Suppose now that det A > 
and that m e . 1 ..' 1 . — .«'c«p+«i = for some sequence e{ , e^ , ••• , e( 2p+ „ . 
Then with d,- = d for all j for which ej = 1 and d,- = d~ l for all j for which 
e; = 0, we have [see equation (18)] det(o7 + J u ) — » as d — » co.n 

3.6 Theorem 3: Let T t 3, /e< F(-) « 3 r o 2p+c> wzYfc eac/i /,-(■) continuously 
differentiable on (—<*>, 00) and /<(a) > /or aZZ a e (— 00, »), a?id feZ 
denote [0, 00) or an interval contained in [0, °o). TVien /or all but at most 
a finite number of points a contained in d, there is a real constant e a > 
such that det((7/ + J u ) ^ e„ for all u e E {2v+q) if and only if det(crl + /„) ^ 
for all a c 4 and aZZ ii e e 12p+ "\ 

3.7 Proof of Theorem 3 

As in the proof of Theorem 1, (err + T 1 )" 1 ^ + (?) e P for all a e 
if and only if det (of + J u ) ^ for all a e # and all n. We shall also use 
the fact that since det (or + T) > for all a ^ 0, each principal minor 
of (or + T)~ l (ac + (7) is a finite- valued rational function of a for all 
a ^ 0. 

(if) If det (a/ + J u ) 9* for all u and all <r s d, then (or + T) '(ac + 
(?) e P for all a 1 H. It is clear that (<rr + T) _1 (<rc + (?) eP for all suffi- 
ciently large <r > 0. Thus each principal minor of (<rr + T)~ (ac + (?) 
is nonnegative for all a t 3 and is positive for all sufficiently large <x > 0. 
They are therefore positive for all but at most a finite number of values 
of a e d. Thus, by Theorem 2, if det(<rZ + Ju) 7* for all a e and all u 
there exist for all but at most a finite number of points rda positive 
constant e„ such that det(<rl + /„) ^ e„ for all u. 

(only if) If det(o-7 + J u ) = for some 0- e £f and some u, then, for 
that a, (or + jT)"'(ac + 6') 4 P . That is, for that a at least one principal 
minor of {or + T)~\<jc + (?) is negative. This means that (or + 7 7 ) -1 
(o-c + (?) ?■ P„ for all o- contained in some interval #' C #, and by Theorem 
2, for all j e d' there is no e„ > such that det(7 + J u ) ^ e„ for all n. □ 
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3.8 Theorem 4: Let R(-) be a continuously differ entiable mapping of 
E" into E", and let J(R) Q denote the Jacobian matrix of R(-) at an arbitrary 
point qt E". If the elements of J(R) Q are bounded on E', and if there exist 
real constants a > and 6^0 such that det J{R) a ^ (a + b \\ q ||) _1 for 
all qe E", then R(-) is a homeomorphism of E" onto E". 

3.9 Proof of Theorem 4 

If Ref. 16 Meyer proves 1 " that R(-) is a homeomorphism of E n onto 
E" if J(R)~* exists for all q e E" and there exist real constants a > 
and > such that || /(fl); 1 || < a + 0\\ q || for all q e FT. 

With q an arbitrary element of E", let Xi , X 2 , ■ • • , X„ denote the 
eigenvalues of J(R)l T J(R) Q , and let X t = min, jX,|. Then \i\ 2 • • • X„ = 
[det J{R) g ] 2 ^ (a + b || q \\)~ 2 , and since the elements of J(R) Q are 
bounded on E", there is a constant X > such that X, ^ X for all j and 
all q t E". Thus 

(X0 1/2 ^ v«'»*-»(a + b\\ q ||)-' (19) 

for all q. For any xzW and any g e E 1 , x tT J(R)l r J(R) a x ^ X^'z; that is, 

|| J(RU || ^ (x.m x || ^ x- a/au - 1} (a + b\\ q ||)-»|| x ||. 
With x = J(R)~ l y in which ?/ is an arbitrary element of E", we have 

\\J(R); l y\\ ^^ 1/2Un ' l) (a + b\\q\\)\\y\\, 

which shows that our hypothesis concerning det J(R) Q ensures that 
Meyer's condition on H/CR),'!! is satisfied. □ 

3.10 Some Further Definitions 

Definition 4: For each T e 3, let 3(T) denote the set of all matrices 71/ 
such that M = M, © M 2 © • • • © M v © /„ with 



M t = 



and 

< 0_ < ,v. 



(*j 



1 


-8l k) ~ 








-—5/ 


1 J 








> and 


8l k> = 


if 


a r = 0, 


> and 


s (t) 

Of 


= 


if 


Ik) n 

«/ = U, 



0^ s (k) ^ (fc) ;r <*) 

< d r ^ a/ it otf 

for all fc = 1, 2, • ■ • , p. As suggested, if q - 0, then M = M, M 2 © 

• • • © M v , while if p = 0, then M = I, . 

f Meyer's result is a generalization of a well-known result of Hadamard. 17 Hada- 
raard proved that R(-) is a homeomorphism if J(R) q ~ l exists for all q t E n and 
satisfies || J{R) q ~ l || ^ a for all q e £" for some positive constant or. 17 
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Definition 5: For each T t 3, let \(T) denote the set of all 2 2p matrices 
M such that M = M x © M 2 © • • • © M„ © /„ with 



and 



1 -|»" 



.-5) 



(*> 



I™ = a™ or Ij" = 0, 

8< k) = a\ k) or 5<" = 0, 

for all k = 1, 2, • • • , ?>. As suggested, if g = 0, then M = M x © M 2 © 
• ■ • © 71/ „ , while if p = 0, then M = /„ . 

Definition 6: Let Q (2pta , denote the family of all 2 (2p+9) - 1 sets w = 
(*i i *a j ■ ■ * >*>■}, including the null set, such that r < (2p + q) and 
wC {1,2, ••• , (2p + g)(. 

Definition 7: For M an arbitrary square matrix of order (2p + q), and 
for each w t Qc 2p+a) , let M„ denote the principal submatrix obtained 
from M by deleting rows t'i , i 2 , ••• , i T and columns i v , i 2 , ■ • ■ , i r • 
(If w is the null set, then M „ = M.) 

Definitions: For each j e (1, 2, • • • , (2p + g)}, let ?7, denote the 
(2p + g) -column- vector with unity in the j'th position and zeros in all 
other positions. 

Definition 9: For each Tt3 and each w t Q (2p+< ,) , let T" denote the 
matrix obtained from T by replacing the jth column of T with [/, for 
all j t w. 

3.11 Theorem 5: Let T z 3. Then the following statements are equivalent. 

(i) M~ X G c P for all M e 3(T). 
(it) (D a + T)~\D b + G) t P for all diagonal D a ^ and all diagonal 

D h ^ 0. 
(Hi) T~\G + D) e P /or aZZ diagonal D ^ 0. 
(«;) (Z)„ + T)" 1 (D b + G) e P /o'" aM diagonal D a > awd a// diagonal 

D b > 0. 
(y) T _1 (G + D) t P /or att diagonal D > 0. 
(vi) (T„) 'G. c P /or aU iu e Q (2p+a) . 
(m) [(r u, ) _1 G] ir c P /or aU u> e Q (2p+0 > . 
(ran) M _1 G e P /or aZ2 flf e 3 n (T). 
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3.12 Proof of Theorem 5 
[(i) and (it) are equivalent] 

By Lemma 1, (D. + T)~\D b + G) e P if and only if det [(D a + T)" 1 
(Z) 6 + G)+D]9*0 for all diagonal Z) >0. Thus (D a + T)~ l (D b + (?) 

e P for all A, ^ and all D b ^ if and only if 

det [(AIT 1 + D a + 7 1 )!) + G] * 

for all D a ^ 0, all D fc ^ 0, and all Z) > 0, and hence if and only if 

det [(A + T)D + G] * 

for all diagonal A ^ and D > 0. Let !T A = (A + T)(I + A)~\ Then 
(£> a _j_ 7 , ) _1 (Z) t + G) e P„ for all D a ^ and all D 5 ^ if and only if 

det[T A (7 + A)D + G] 9* 

for all A ^ and all D > 0, and hence if and only if det (T K t> + G) ^ 
for all diagonal Z) > and all A ^ 0. By Lemma 1, this means that 
T~ K 'G i P for all A ^ if and only if (D a + T)'\D b + G)zP for all 
D a ^ and all D„ ^ 0. We observe that T A = (T A ) t © (T x ) 2 © ■ • • 
(^a) p © I« in which, with A = diag (X, , X 2 , • • • , X (2p+o) ), 



(T K ) k = 



1 



(A) 



(A.-) 



.1 + X 2ft -i 



1 + x 2 , 

1 



for fc = 1, 2, • • • , p. Thus for each A ^ 0, T A t 3(5P); and if il/ is an 
arbitrary element of 3{T), there is a A ^ such that M = T A . There- 
fore (D a + rr'CD* + G) e P for all D a ^ and all Z) 6 ^ if and only 
if M _1 G e P for all M z 3(T). 

[(i) and (m) are equivalent] 

Repeat the proof of "(i) is equivalent to (ii)" with each statement 
that Z) a ^ replaced with D a - diag (0, 0, • • • ,0). 

[(ii) and (/y) are equivalent and (Hi) and (y) are equivalent] 

Suppose that (ii) and (iv) are not equivalent. Then (D a + T)~ l (D b + 
G) c P for all D a > and all Z) 6 > 0, and for some D* ^ and some 
-D* fc ^ 0, with D* > 0orD% > or D* > and D* b > 0, (D* + 
T)-\D% + G) * P . Thus some principal minor of (D* + ZT^t + G), 
and hence of (D* + T)~ l (D* h + G) det (D* + T), is negative. Let 
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m(D* , D%) be some negative principal minor of (D* + T)~ 1 (D% + 
G) det (D* + T), and let m(D* + el, D* b + el) be the corresponding- 
principal minor of (D* + el + 7 , )~ 1 (D* 6 + eJ + G) det (D* + e/ + T) 
for all real e ^ 0. Thus m(D* + e/, D* b + ei) is a polynomial p(e) in e 
for e ^ 0, and p(e) ^ for all e > 0. Therefore p(0) ^ 0, which con- 
tradicts m(D* , D%) < 0. 

A proof that (Hi) and (v) are equivalent can be obtained by modifying 
the previous paragraph in an obvious manner. 

[(vi) is equivalent to (v)] 

By Lemma 1, T~ l (G + D) tP, for all diagonal D > if and only 
if det \T-\G + D) + £>*] ^ for all diagonal D* > and D > 0, and 
hence if and only if det (G + TD* + D) ^ for all D* > and all 
D > 0. Therefore, by Lemma 1, T _1 (G' + D) c P for all i) > if and 
only if (G + TD*) z P„ for all D* > 0, that is, if and only if det [G w + 
(TD*)„] ^ for all w z Q {2v+a) and all D* > 0. Since (TD*) W = T U D* , 
we see that T~\G + D) z P for all D > if and only if 

det[(T w )~ 1 G w +Dt] ^0 for all u'cQ (2p+a > and all Z>* > 0. (20) 

But, by Lemma 2 (which follows) condition (20) is equivalent to the 
condition that det [(T w y 1 G w + Z>*] > 0, and hence that det [(T a y'G w 
+ D*] t* 0, for all w z Q (2p+a) and all D* > 0. Thus by Lemma 1, T'\G 
+ D) z P for all D > if and only if (T, i y l G w z P for all w z Q (2p+s) . 

Lemma 2: If A is a real square matrix of order n such that det (D + A) = 
/or some diagonal D > 0, Men det (/) + ^4.) < /o?- some diagonal D > 0. 

Proof: Using the notation of the proof of Theorem 2, 

det (D+ A) = £' </;• £' • ■ <C m*,.......... (21) 

for all Z) > 0. Since m ,.,....., = 1, if det (D + .4) = for some D > 0, 
then for at least one sequence e[ , e', , ■ ■ ■ , e', we have m,., .,'.,.....», < 0. 
If ?n .o..--.o = det i4 < 0, then there exists a positive constant <x, such 
that det (D + A) < whenever < c/, < a^ for all ;'. If det A ^ 0, 
then, with d, = d for all ; such that e] = 1 and d, = d _1 for all ; such 
that e'f — 0, there exists a positive constant <r 2 such that det (D + A) < 
for all d > a 2 [see (21)]. □ 

[(trc) and (y») are equivalent] 

We shall prove that 

[{TT'G]* = {T w r l <~;«. for all w e Q (2p+a) . (22) 

Obviously the equality of (22) is satisfied if mj is the null set. 
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It is convenient to introduce the following notation. Let u denote the 
1X1 matrix containing the entry 1. Let <p denote what might be called 
the empty matrix, a matrix with no rows or columns; by this we mean 
that <p is to be interpreted in the following manner: <p <p — <p, I, — <p 
when s = 0, (p' 1 = (p, and if M x and M 2 are any two (ordinary) matrices, 
then <p © M t = M x , M x © <p = M , , and M x ® ? © M t = ¥,© M a . 

Let w e Q( 2 p+,) and let w not be the null set. The matrix T can be writ- 
ten as the direct sum 7\ T 2 • • • T v © /, . In terms of w and ip, 
T a = 1 1 © t 2 © • • • © t p © I, , in which s = q — q where g is the num- 
ber of elements contained in the intersection of the sets w and ( 2p + 1 , 
2p + 2, • • • , 2p + g] , and for k = 1, 2, • • • , p: t k = T k if both (2k - 1) 
and 2k are not elements of w, t k = <p if both (2fc — 1) and 2k are elements 
of to, and t k = u if either (2fc — 1) e to and 2/c e" w or (2fc — 1) «" to and 
2k e to. Thus (r^)" 1 = C © £* © • • • © C © /. . But (TT 1 = 
^r 1 ©r 2 - 1 ©---©f'; 1 ©/ a ,inwhichforA : = 1,2, • • , p: t „ -- T k 
if both (2k — 1) and 2k are not elements of w, 

"l 
1 

if both (2k — 1) and 2k are contained in id, 

1 



rr - 

if (2/c — 1) e w and 2fc e" w, and 

t: 1 = 










\-CXf 



1 



if (2k - 1) I w and 2fc e m>. Thus we see that [(T*)~ l ] m = (T.)" 1 . Let 
(U , ) (7 7U ') -1 denote the (2p + g — r) X (2p + g) matrix obtained from 
(T*)' 1 by deleting rows ii , i 2 , ■ • ■ , i* • But all elements of columns 
*i > *2 j ■ ■ • , ir of ( „)(jP") -1 are zeros, and hence, with (7(„) the matrix 
obtained from G by deleting columns i t , i 2 , ■ ■ ■ , « r , 

[(wm) and (i) are equivalent] 

If M~ l G t P for all M t 3„(7 1 ), then [(T w y l G] w t P for all w z Q (2p+Q) . 
Thus, statement (viii) implies statement (vii). Since we have proved 
that (vii) is equivalent to (i), it suffices to prove that (i) implies (viii). 
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Suppose that M~ l G t P n for till .1/ e 3(T). Let M be an arbitrary ele- 
ment of 5 (T). Then [M + S(T - M)] t 3(T) for all 5 t (0, 1], and 
therefore [M + S(T - M)\~ l G t P for all 8 z (0, 1]. At this point a 
continuity-type argument similar to that used in the proof of [(it) and 
(iv) are equivalent] shows that M~ l G e P . □ 

3.13 Corollary 1 (Corollary to Theorem 5): 

IfTtZ and M~ l G t P jor all M t 3(T), f/ie?i det (<rl + J u ) 5* /or aZZ 

<r ^ and oM it c E (2v+Q) provided that for all j /,(•) is continuously dif- 
ferentiable on(-», <*>) and /J (a) > for all a e (— co , co). 

3.14 Proof of Corollary 1. 

If T e 3 and M _I G e P„ for all Af e 3(T), then, by the equivalence of 
(i) and (w) of Theorem 5, (ct + T) _1 (c7C + G) c P for aU a ^ 0. The 
first portion of the proof of Theorem 1 shows that if (ar + T) _1 (crc + 
(?) e P for all a ^ and if for all j /,(•) is continuously differentiable 
on (— oo , cc ) and f'(a) > for all a e (— °o , °o ), then det (<rl + ./„) ^ 
for all a ^ and all u t E V2p * Q) . 

3.15 Definition 10: For p > let Q( 2p+ „) denote the subset of Q (2 p +0 ) 
containing all sets 10 belonging to Q( 2p+a) such that w is not the null set 
and 2k e w if and only if (2fc — 1) e if for fc = 1,2, • • • , p. For p = 0, 
let Q[ 2v+a) denote the family of all sets contained in Q (2p+9 ) with the 
exception of the null set. 

3.16 Theorem 6: If T e 3 and T~ l G t P , then (T„)~ l G a t P for all 
w t Q[ 2v+Q) ■ 

3.17 Proof of Theorem 6 

Let T e 3, and let T~ l G e P . By Lemma 1, det (TD + G) ^ (and 
hence det (TD + G) > 0) for all diagonal D > 0. Let w = {i x , i 2 , • • • , 
*rl e Q(2p+,) » a "d let d, t = cZ for ft = 1, 2, • • • , r. 

It may be the case that (TD + G) is a block matrix of the form 



H 2] («T + tf 22 )_ 



(23) 



in which f is a direct sum of all 2 X 2 and 1X1 block matrices on the 
diagonal of T which do not appear in T w , and H 12 , H 21 , and H 22 are 
independent of D. Clearly det f > 0. If (TD + G) is not of the form (23), 
then by a sequence of interchanges of rows and corresponding columns of 
(TD + G) we obtain a matrix of that form. 
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Thus, for some t of the form indicated above and for the correspond- 
ing constant matrices H 12 , H 21 , and H 22 whose elements are elements 
of G, 

det(Tn + G)=det[( TD +V» (c /+tfJ 

for all cZ, > for j 4 w. For all sufficiently large d > 0, det (df + H 22 ) 
> 0, and then 

< det (TD + G) = det (df + tf 22 )-det [(TD + G) w 

-H 12 (df + H 22 y 1 H 21 ] 

for all dj > for j 4 w. Since H l2 (dT + H 22 )~ l H 2 i approaches the zero 
matrix of order (2p + q — r) as d —> °° , we must have det (TD + G) w 
^ for all di > for j 4 w. Therefore, since (TD) W = T W D W , we must 
have det (T W D W + G w ) ^ for all D w > 0. But this means (see Lemma 
2) that det (T W D W + G w ) ^ for all D w > 0. Thus, by Lemma 1, 
(T w y l G w zP . □ 

3.18 Theorem 7: If T e 3 iw'tfi p = 1 or p = 2, and t"/ T~ l G e P witfi 
(7 the short-circuit conductance matrix of a transformerless ■positive- 
element resistance network, then (T w )~ l G w e P for all w e Q (2p + ) • 

3.19 Proof of Theorem 7 

Suppose that T' X G e P with p = 2. Theorem 6 asserts that (T„)~ l G w 
e P Q for all w e Q[ 2p+Q) . But, aside from the null set, the sets 
w = |t'i , i 2 , • ■ ■ , i T \ that are contained in Q<, 2p + Q ) but not in Q{ 2pi . Q) 
possess the property that T w — Ti ® I i2 + Q - r) , or T w = u 7^ 
/ (1+(( _ r) where ij is the 1X1 matrix containing the element 1, or 
T = I 

If T w = J (4+a - r) , then obviously (T w y 1 G w z P . U T w = T,0 
/ ( 2 +fl - r ) , then for any D w = diag [D 2 D (2+a _ r) ] with D 2 > and 
Z) (2+a _ r) > diagonal matrices of order 2 and (2 + g — r) respectively, 



det (T W D W + <?„) = 



T,D 2 + G ll (? 12 

G 2i Z) (2+a _ r) + (7 22> 



(24) 



in which G n , G l2 , G 2i , and G 22 are the appropriate block matrices 
of G w . Since det [Z) (2+a _ r) + G 22 ] > 0, we have 

det (T W D W + G w ) = det [Aa+,-r) + G' 22 ]-det [T 1 D a + G'„ 

— Gi 2 [i)( 2 4-a-r) + G22] G' 2 ,j. 
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But G u — G 12 [D i2+a -r) + G 2 2]~ 1 G 2 i is the short-circuit conductance 
matrix of a transformerless common-ground 2-port network; it is of the 
form 

_ — 012 022 J 

with U ^ 0, g 22 ^ 0, </, 2 ^ 0, r/ n ^ </, 2 , and g 22 ^ 12 . Therefore 1 

det {T 1 D 2 + (?„ - (?i 2 [A 2+0 -rt + G 22 ] _1 (7 21 | > 

for all D 2 > and all D (2+0 _ r , > 0, det (T W D W + G w ) ^ for all £„, > 0, 
and hence, by Lemma 1, (T,,,) -1 ^ c P . Finally, the case in which 
T w = u T 2 /(i+ a -r) can be treated in a manner similar to that used 
to show that (T w )" l G w £ P when T,„ = 7\ / (2+< ,- r) , since, with w 
such that T w = w T 2 7 (1+(2 _ r) , and with D an arbitrary diagonal 
matrix of order (4 + q — r), a sequence of interchanges of rows and 
corresponding columns of (T W D + G w ) can be performed to obtain a 
matrix of the type that appears on the right side of equation (24). 
Therefore {T w y x G w e P for all w e Q (2p+a) . 

When p = 1, aside from the null set, the sets w = [ii , i 2 , • ■ • , i r ) 
that are contained in Q i2p+a) but not in Q[ 2p+Q) possess the property 
that T w = / (2 ^_ r) and obviously when T w = J (2 + a - r , , (T 7 ^)" 1 ^ eP . □ 

3.20 Theorem 8: Let T e 3 and let G possess the property that -for some 
diagonal matrix D > 0, both DT and DG are strongly-column-sum dom- 
inant. For each j = 1, 2, • • • , (2p + Q) let /,(•) be a continuous mono- 
tone-nondecr -easing mapping of E l into itself such that /,(0) = 0, let 
h t (0, co), and, with F(-) and C(-) defined relative to the /,(•) as in 
Section 2.1, suppose that the sequences {y n \ and {w n } in E l2p+Q) satisfy 

y n+1 + h{TF[C-\y n+1 )] + CC" 1 ^,)} = y n + w n 

for all n ^ 0. Then there exists a positive constant 5 depending only on the 
c, , the t,- , T, G, and D such that 

(i) || Dy n ||, S (1 + 8h)- || Dy ||, + £ (1 + SA)'* || Dw {n - k) ||, 

*=i 

/o?- all n ^ 1 , and 

(ii) 1 1 D(t/„ - Vn) I Ii ^ (1 + 5/i) ""I | Z%o - £o) 1 1. +£ Z(l + 5/i)~ fc 

k = 

for alln ^ 1, in which \y n \ is any sequence in E i2v * a) with the property 
that \\ D(y„ — y* n ) || x ^ e for all n ^ 1 with e a positive constant and the 
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sequence |?/*| such that 

y* +1 + h[TF[C-\y* + i)] + GCT'ty&J} = y„ + w„ 
for all n ^ 0. 

3.21 Proof of Theorem 8 

We shall first prove part (n). With D such that DT and DG are 
strongly-column-sum dominant, we have for all n ^ 

Dy n+1 + h\DTF[C-\y n ^)] + DOtT'^OI = Dy n + ZM, 

and 

%„* + i + h{DTF[C-\y* + i)] + DGC~\y* +1 )} = Dy* + JD(fif. - */*) + Dw n 

in which we shall take y% to be y . As in the proof of Theorem 2 of 
Ref. 3, we write 

F[C-\y n+1 )] - F[C-\y* +l )\ = diag ( - + K ^ )y ) to*. " I&0 (25) 
and 

C~\y n+1 ) - C'\y* +i ) = diag ( - - ^J fo+i ~ tfW (26) 

in which r(w), depends on the jth components of y„ + , and t/* +1 , and 
r(n)j ^ for all n ^ and all j. 

Thus, with Q = DTD' 1 and fl = DGD~\ 

{' + ^ diag fcrfe) + w diag t+^)K- ~ *" 

= Dfa. - ?/*) - £>(£„ - yt) 

for all n ^ 0. At this point we shall use the proposition that if M is any 
real matrix of order (2p + q) with the property that there exists a 
positive constant rj such that m /y — ^ lV , | m,,- | ^ 77 for all j, then 
|[ Mx ||t ^ 7j || a ||x for aU x e # l2 " +a> . Now let 

M = ll + hQ diag ( , r(n) ' / r ) + W2 diag ( , * , > ) 

i & \c,- + r,r(7i),7 & \c ; + T,T{n)il 

for arbitrary n ^ 0. Then for arbitrary j 

.. - £ , m ,.., = 1 + fe gi ./ r(w) ' ) + hrJ l ) 

f£ \Cj + Tjrinjj/ \Ci + T,r(n)j/ 

V"^ r(tt)f 1 

— h ?. Qij : ;r— + Ta j TV 



m 
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in which 

5 = min S mm e^ta - 2] | r, y | I , min rjHqu - 2 I Qa 

Therefore 

II #(«/..+ 1 - 2/?+0 Hi 

< (1 + 5/t)" 1 || D(y„ - y*) - D(y n - y*) ||, 

S (1 + a*)" 1 || D(2/» - ?/*,) Id + (1+ 5/ir 1 || D(j?„ - 7/t) II. 
g (1 + 5/i)" 1 || Dfo„ - y*) |„ + e(l + ShT 1 
for all ?i ^ 0, and hence 

|| D(y n - y*) ||, * (1 + 5/i)- || Z)(y„ - i/o*) ||i + eE(l+ SZ0~* 

A-l 

for all n ^ 1. Finally, since \\ D(y n - y n ) ||, ^ || Z%,. - 2/*) II, + 
I! D(y* — y n ) ||i ^ || Z%„ — y* n ) ||i + e, and since y% = y , 

II D(y n - y n ) ||i £ (1 + «*)"" II D Ga> - ft) ||i + «S (1 + Sh)'" 

for all n ^ 1, which completes the proof of part (ii) of the theorem. 
The proof of part (i) is similar to that of part (ii) . Using 

Dy n+l + h{DTF[C-\y n ^)] + DGC -1 ^,)) = -%„ + JDtft. 

for all ?i ^ 0, and equations (25) and (26) with y* +1 = 6 for all n, we 
find that 

|| Dy n+1 |K g (1 + ft 5)" 1 II £h/n ||: + (1 + h 8)- 1 || D«>„ ||, 

for all n ^ 0. Therefore 

II Dy„ ||i ^ (1 + &«)"" || Z>2/ ||, + S (1 + h8)- k || Dw (n - k) ||. 

it-i 

for all n ^ 1. D 
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3.22 Theorem 9: Let T z 3 and let G possess the property that for some 
diagonal matrix D > 0, both DT and DG are strongly-column-sum dom- 
inant. Let B(-) denote a real continuouslij-differentiable (2p + q)-vector- 
valued junction of t for 1 1 [0, °°) such that both B(-) and (d/dt)B(-) are 
bounded on [0, «>). With F(-) such that each /,(0) = 0, and ivith C(-) 
defined relative to F(-) as in Section 2.1, let u(-) satisfy 

%} + TF[C~\v)] + GC'\u) = B(t), t ^ 
at 

and, with h an arbitrary positive constant, let u„ denote u(nh) for all 
n ^ 0. Let \y n \ be a sequence in E l2p+Q) such that 

y n+1 + h\TF[C^(y a+1 )] + GC~\y n+1 )] = y n + hB[(n + l)h], n ^ 0. 

Then there exist positive constants 8 and p, both independent of h, such that 

II D(u n - y n ) ||x £ (1 + 8hy n || D(u - y ) ||i + P h 

for all n ^ 1. 

3.23 Proof of Theorem 9 

The sequence {u n } satisfies 

u n+1 + h{TF[C~\u n+l )] + GC-\u n+1 )} 

= u n + B[(n + l)h] + kn , n^O 

in which £„ is often referred to as "the local-truncation error at step 
n." We shall first bound £„ . 

Since B(-) is bounded on [0, <*>), and since for some D > 0, both 
DT and DG are strongly-column-sum dominant, a direct modification 
of the proof of Theorem 1 of Ref. 5 shows that u(-) is bounded on 
[0, oo ); and hence since 

~ = J u {TF[C-\u)\ + GC-\u)} - JJS{t) + J- 5(0, / ^ (27) 

with (d/dt)B(-) and the elements of the Jacobian matrix J u bounded, 
it is clear that (d 2 u/dt 2 ) is bounded on [0, «>). By the usual Taylor- 
series-type argument we can show that for arbitrary n ^ 0, £„ = %h 2 U n 
in which for each ; the jth component of U n is the jth. component of 
{d 2 u/dt 2 ) evaluated at some point contained in the interval [nh, (n + 
l)h]. Thus there exists a positive constant p v such that 

||%||i ^ ¥i 2 Pi for all n ^ 0. (28) 
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Therefore, using (28) and the equations 

m, 1+1 + MTW-V+,)] + GC'\u„^)\ 

= u n + B[(n + l)h] + £„ , n^O 

Vn+i +h\TF[C- 1 (y n > l )} + GC- J (y n+i )\ 

= y n + B[(n + 1)*], n ^ 

by an argument similar to that used in the proof of part (ii) of Theorem 
8, and with 5 as defined there, we find that 

II *>(«.+, - 2/n + II. ^ (1 + ft)"' || D(u n - y n ) ||, + (1 + ft)- 1 }*'* 
for all n ^ 0, and hence that 

1 1 D(u n - y n ) ||, g (1 + 5/i)-" || D(u - y ) \\ t + jA 2 Pl E (1 + *&)"* 

£ (1 + «ft)~" || Z)( Mo - 2/o) 111 + Ife'p. Z (1 + «fc) _ * 

fc-l 

^ (1 + 5/0~" || D(w - y„) II. + hhS~ 1 Pl 
for all n fc 1. □ 

3.24 Definition 11: Let It = #, © tf 2 © ■■■ ® R P © R in which 
fl = diag ixx , r 2 , • • • , r a ) with ?•, ^ for j = 1, 2, • ■ • , q and 



Rk = 






with rj w ^ 0, r l b k) ^ 0, and rj k) ^ for all fc = 1, 2, • • • , p. As suggested, 
if q = 0, then # = fl, © fl 2 ©•••©#, , while if p = 0, then R = R . 

3.25 Theorem 10: Let T e 3. If p > and if R satisfies 

aj w (l - 0"*V. (W = rj 4) 

«J M <1 - a l f k) rVc k} = r ( b k) 

for k = 1, 2, • • • , p, then T~ X G{1 + RG)~ X e P whenever T~ X G t P . 

3.26 Proof of Theorem 10 

By Lemma 1, T l G(l. + AG)" 1 e P if and only if 

det [T~ l G(I + AG)" 1 + D*\ * (29) 

for all diagonal D* > 0. But (29) is satisfied if and only if 
det (T~'G + D*RG + Z)*) * 0. 
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Here, since 

a ( r k) (l - O'VJ* - rl k) 
«}»(! - a ( f k) yV c k} = r[ k) 

for k = 1, 2, • • • , p we have # = DT -1 for some diagonal matrix D ^ 0. 
Thus (29) is satisfied if and only if 

det [(/ + DD*)T~ l G + D*\ * 0. 

When T X Q z P we have 

det (T' l G + D) 5* 

for all diagonal D > 0. Thus (29) is satisfied for all D* > whenever 

T-'GtPo. n 

3.27 Theorem 11: If M~ X G z P for all M z 3, then for any T z 3, 
Z 1-1 ^/ + AG')"" 1 £ P for all R. 

3.28 Proo/ of Theorem 11 

Let T e 3. As in the proof of Theorem 10, T~ l G(I + RG)~ l z P if 
and only if 

det [(T- 1 + D*fl)G + Z)*] ^ 

for all diagonal D* > 0. It is a simple matter to verify that for each 
D* > and each R there exists an M z 3 and a diagonal matrix D > 
such that (T -1 + D*R) = DM' 1 . Since iV/ _1 G z P for all M c 3, we 
have (by Lemma 1) 

det (Dfr'G + D*) * 
for all D* > 0. □ 

3.29 Theorem 12: Let T z 3 wi% p > and q ^ 0. 27ien M~ l G c P„ /or 
aZ/ M z 3(T) # and o?% i/ T _1 G(7 + PG) _1 z P /or a22 # such that 

al k \l - 0"Vi w ^ r* M 

a < fc) (i - o-v*" ^ ri M 

/or fc = 1,2, ■ • • , p and r,- ^ /or a22 j such that 1 ^ j ^ o. 

3.30 Proo/ o/ Theorem 12 

As in the proof of Theorem 10, T' l G(l + RG)~ l z P„ if and only if 
det (T~ l G + D*RG + D*) ^ (30) 
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for all diagonal D* > 0. The inequalities r, ^ for all j such that 1 ^ 
j ^ q and 

a r (1 — a T ) r e ^ r b 

a/ (1 — a/ J r e ^ r 6 

forfc — 1, 2, • • • , p are equivalent to the condition that R = D v T~ l + D 2 
for some diagonal matrix Z> 2 = and some diagonal matrix D t t S, in 
which £ is the set of all diagonal matrices D ^ such that DT' 1 is 
symmetric. Hence T~ X G{I + RG)~ l t P for all such i? if and only if 

det ![(/ + DJ)*)T X + D*D 2 ]G + D*\ ^ (31) 

for all diagonal D* > 0, D 2 ^ 0, and Z), e 5. 

Let A = diag (X, , X 2 , • • • , X (2p+0) ) be such that 

D 2 = D*- 1 A _1 A(7 + D X D*) 

in which 

A = diag (5, , 5, , o 2 , o 2 , • • • , o„ , 8„) 7, 

if g > 0, A = diag (5, , 8, , 8 2 , 8 2 , • • • , 8 P , 8 P ) if q = 0, and 

5, = 1 - aJV w for k = 1, 2, ■ ■ • , p. 

The left side of (31) is 

det [(/ + i) 1 Z)*)(7 , " 1 + A~\\)G + D*] 

which can be written as 

det [(/ + DJ)*) A~\I + A) A,T~ K l G + D*\ (32) 

with 



and 



7V = AX 1 A(7 + A)- l (r _I + A-A) 



A A = diag (8[, 81, 8',, 8',, ■■■ , 8 P , 8' v ) © L q 



if q > and A A = diag (6( , 5( , 5£ , 8' 2 , ■ ■ ■ , 8' p , 8' p ) if q = 0, in which 
forfc - 1, 2, ••• , v 

Si = 1 - a}V M (l + Xc,*-!,)- 1 ^ + X,,)" 1 . 

But (32) vanishes if and only if det (T~ K l G + D) vanishes, in which 
D = A"'(7 + A)" 1 A(7 + DiD*)'^*. We observe that D is a positive 
diagonal matrix and that given any diagonal D' > and given any 
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A ^ we can choose D* > and D, e S so that D = D'. Thus T l G(I + 
RG)' 1 t P for all R = (D 1 T~ 1 + D 2 ) with D, e S and Z> 2 ^ if and 
only if 

det (Tl l G + 5) 5* 

for all A ^ and 7) > 0, that is, if and only if T?G t P for all A ^ 
(see Lemma 1 of Section 3.1). But 

T A = Tx © T 9 © ■ • • © T, @ I, if q > 



and 
with 



T A = 7\ © T 2 © • ■ • © T p if g = 



T k = 






(*) 

— Or 

1 + X»_, 

1 



_1 + x 2J 

for all fc = 1, 2, • • • , p. Therefore T' l G(I + RGy 1 e P for all # = 
(DjT -1 + D 2 ) with D 2 ^ and A e S if and only if M l G z P for all 
Mt 3(7). □ 

3.31 Definition 12: Let J3 denote the set of all F(-) such that 

(t) nOeffi'^'.and 
(m) for each j = 1, 2, • • ■ , (2p + g) there exists a real constant 

/3, such that /,(•) is a strictly-monotone-increasing mapping of 

E 1 onto either (j8 y , «>) or (— «, 0,), and 
(m) whenever p > 0, /<»-«(•) and /»(•) are both bounded on either 

[0, 00) or (— * , 0] for k = 1,2, •■■ , p. 

3.32 Theorem 18: Let T c 3, and, referring to the network of Fig. 1 m 
which it is assumed that R (see Section 2.1) is the zero matrix, let G denote 
the short-circuit conductance matrix of the linear portion of the network. 
(The linear portion is assumed to contain only sources and linear resistors 
of nonnegative resistance.) Then the equation F(x) + T~ l Gx = B pos- 
sesses a unique solution x for each F(-) £fF 3 and each B zE " +0 if and 
only if T~'G t P Q and det G 9^ 0. If T~'G c P and det G = 0, then there 
exists a real (2p + q)-vector jj such that (i) 77 5^ 0, and for some F(-) c $3 
all of the components of F(a-q) are bounded on a e [0, »), and (ii) for any 
F(-)z5 a with the property that all of the components of F(ar)) are bounded 
on a t [0, 00) the equation F(x) + T~ x Gx = B does not possess a solution 

/Ti 77»(2p + g) 

or some B e 1L 
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TRANSISTOR 1 



TRANSISTOR p 



L (2p-l) 



-ap 



+ v,- 6 -v 2 + 




o +v ( 2p -i)-6 ~ v 2p + 6 



L (2pH 



*C2p+1) 



L (ap + c^) 



v t2p + q)" 



diode q 



Fig. 1 — General network containing transistors, diodes, resistors, and sources. 

3.33 Proof of Theorem 13 

(if) If T~'G t P with det T~ l G t± 0, and if ¥{■) e SF a , then, since 
each /,-(•) is a strictly-monotone-increasing mapping of E l onto (j8,- , °°) 
or (— oo, /3,-) for some real constant /3, , by Theorem 4 of Ref. 2, the 
equation F(a:) + T~ i Gx = B possesses a unique solution x for each 
BtE (2p+,1) . 

(only if) Assume that T~ l G 4 P . Then since fF 3 is contained in 
$o 2p+< " » by Theorem 1 of Ref. 3, for each F(-)z5 3 there exists aBe E i2v+Q) 
such that there are at least two solutions x of F(x) + T~ x Gx = 5. 

Assume now that T^'G t P and that det G = 0. We shall use the 
proposition that if R(-) is any continuous mapping of E i2p+a) into itself, 
then R(-) is a homeomorphism of E i2p+a) onto itself if and only if R(-) 
is a local homeomorphism on E l2p * Q) and || R(x) || — > « as || * || — » <».t 

Let fl(-) be denned by the condition that fi(x) = F(:c) + r _1 (?x for 
all x t E Vip *" ) . For any F(-) t JFg the operator #(•) is a local homeo- 
morphism on ^ <2 " +0) f since with F(-) such that each /,(•) is a strictly- 
monotone-increasing mapping of E x onto E l the mapping [F(-) + T~ G\ 
is a homeomorphism of E Vip * Q) onto itself. 1 In addition, for any F(-) t IF3 
and any B e E V2p+Q) , there is at most one x t # )2p+0) such that R(x) = B. 1 

Let us suppose that for each B t J5 <2p+a) and each F(-) t {F 3 there 
exists a solution x of R(x) = B. Then for all F(-) t £F 3 , R(-) is a homeo- 
morphism of E l2p + Q) onto itself, and hence for all F(-) t JF 3 || R(x) || -> 00 
as || £ || — > °°. But, by Lemma 3 (which appears below) E (2p+Q) contains 
a vector t? such that rj 5* d, jj,- c {0, +1, — 1] for all ;', and Gtj = 0; and if 



t See Ref. 12 and the appendix of Ref. 13. 



1774 THE BELL SYSTEM TECHNICAL JOURNAL, OCTOBER 1970 

p > 0, 7] satisfies rj (2 k-i)V2k ^ for all lc = 1, 2, • • • , p. Let 5 3 (y) denote 
the subset of SF 3 containing all elements F(-) with the property that 
fiiarif) is bounded on a t [0, ») for all ; = 1, 2, ••• , (2p + ?)• Since 
r)(2k-i)V2k ^ for all fe = 1, 2, ■ ■ ■ , p when ;> > 0, it is clear that ^ 3 (v) 
is not empty. However, for any F(-) t JF 3 (i7) we have || R(ari) || = 
|| F(arj) || with || F(ai]) |j bounded on a t [0, <»), which contradicts the 
assumption that there exists a solution x of R(x) = B for each F{- ) t fF 3 
and each B e £ (2p+0) . 

Lemma 3: Let G be the short-circuit conductance matrix of the linear 
portion of the network of Fig. 1. // det G = 0, then there exists a vector 
v e E i2p+q) such that G v = 6, r, ^ 6, and t?, e |0, +1, -1( for all j = 
1, 2, • • • , (2p + q); and if p > rj also satisfies y l2 k-»V2k ^ for k = 
1, 2, • • • , p. 

Proof of Lemma 8: 

Let N denote the (2p + ?)-port resistor network obtained from the 
network of Fig. 1 by removing all transistors and diodes and by setting 
the value of each source to zero. The short-circuit conductance matrix 
G possesses the property that if v e j? (2p+a) denotes the vector of port 
voltages of iV and i e E l2p+Q) denotes the corresponding vector of port 
currents (with polarities as indicated in Fig. 1), then i = —Gv. 

Let det G = 0. Then the open-circuit resistance matrix of N does not 
exist. Therefore there exists a port I of N such that there is no path 
through resistors of N that connects the two terminals of port I when 
all other ports are open-circuited. Let a one- volt source be placed at 
port I so that v ( = 1. Then when all ports j of N with j ^ t are open- 
circuited, i t = and there is zero current in every resistor of N. Let S 
denote a set of port numbers of N with the following properties. The 
number I is not contained in S and when all ports j with j t S are short- 
circuited and all ports j with j $ $> VJ j^| are open-circuited then zero 
current flows through the one- volt source; when any port ;\ f! S W l^l 
and all ports j with j £ S are short-circuited and all ports ; with j «" S U 
\l, ji] are open-circuited then nonzero current flows through the one- 
volt source. It is clear that such a set S exists (with the understanding 
that S might be the null set). In general S contains r port numbers 
where ^ r ^ (2p + q — 1). 

If r = (2p + q — 1), then with v ( = 1 and with all remaining com- 
ponents of v equal to zero, we have Gv = 0. Obviously in this case we 
can take the vector -q of the statement of Lemma 3 to be v. 

If r ?* (2p + q — 1), then, with v, = 1, with v t = for all j c S, 
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and with all ports j $ S \J \(\ open-circuited, there exists for each 
j i S W \t\ some path through the one- volt source and the resistors of 
N that connects the two terminals of port j. Therefore when r ^ (2p + 
q — 1), when all ports j 4 SU \l\ are open circuited, when v t = 1, and 
when V; = for all j e S, the open-circuit voltage »,- at each port j with 
j $ S VJ \t\ is well defined and nonzero. Since no current flows in any 
resistor of N when v t = 1, i>, = for all j e S, and all ports j i SVJ \l\ 
are open-circuited, it follows that v,- e |— 1, +1} for all ; i S. With 
v, = 1, with v,- = for all ;* e S, and with i>, the corresponding open- 
circuit voltage for each j i S U (/j,we have Gv = 9. When p > 0, the 
vector v also satisfies the condition that i^i-i)^* ^ for all k = 
1, 2, • • ■ , p since if ^m-n^t- were negative for some k, then for that 
k y (2 A-i) = 1 and r 2t . = — 1 or i\ 2k _ t ) = — 1 and ?% = 1; in either case 
| f<a*-i) — ^2* | = 2 which contradicts the proposition that a network 
of nonnegative resistors can have no voltage gain. □ 

APPENDIX* 

A theorem due to R. S. Palais 1 " asserts that if R(-) is a continuously- 
differentiable mapping of E" into itself with values R(q) for q e E" , 
then R(-) is a diffeomorphism 1 of E n onto itself if and only if 

(i) det J„ 5* for all q e E n , in which J a is the Jacobian matrix of 

R(-) with respect to (/, and 
(ii) || J?(g) || — > «5 as || g || — > «=. 

If R(-) is any twice-continuously-differentiable mapping of E' into 
itself such that conditions (i) and (ii) of Palais' theorem are satisfied, 
then E n contains a unique element x such that R(x) = 8 in which 6 
is the zero element of E 1 , and there are steepest decent as well as 
Newton-type algorithms each of which generates a sequence in E n that 
converges to x. To show this, let 18 f(y) = || R(ij) || 2 for all y e E n in 
which ||-|| denotes the usual Euclidean norm (i.e., the square-root 
of the sum of squares). Since condition (J) of Palais' theorem is satisfied, 
the gradient V/ of /(•) satisfies (V/)(</) s* 6 unless f(y) = 0, § and 
since condition (ii) of Palais' theorem is satisfied, the set S = \y e JET : 
f(y) ^ f(x i0) ) } is bounded for any x m t E n . Therefore we may appeal to, 
for example, the theorem of page 43 of Ref. 18 according to which for 
any x l0) t E n , for any member of a certain class of mappings <p(-) of S 

* The material of this appendix together with some misprints appears in Ref. 3. 
t See Ref. 12 and the appendix of Ref. 13. 

* A diffeomorphism of E n onto itself is a continuously differentiable mapping of 
E n into E n which possesses a continuously differentiable inverse. 

* Here we have used the fact that ( VfK'j) = 2J v tr R{y) for all y e E n . ls 
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into E", and for suitably chosen constants y„ i 7i , • • • j the sequence 

x m , x (1) , ■ ■ • denned by 

r (. +1) = x <« + 7rf ,( a .c» ) for all k ^ 

belongs to S and is such that \\ R(x (k ') || — > as A; — > «. However, 
since -R _1 (-) exists and is continuous, it follows from 
x- (fc) = R- l [R(p w )] for all fc £ 
and the fact that R(x lk) ) — > as fc — > oo, that lim i _ og x <w exists and 

limas™ = R'\d), 

k—OB 

which means that lim A ._ ro x U) is the unique solution x of R(y) = 0. 
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